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ABSTP.ACT 

Bifurcation theory ia used to analyse the nonlinear dynamic stability 
characteristics of ,:n aircraft subject to single-degree-of-frecdem , 0.3., 
pitching-or rolling-motion perturbations about its trussed steady flight. 

The requisite moment of the aerodynamic forces in the equations of notion 
is shown to be represent rblc in a fora equivalent to the response to finite- 
amplitude osv'illations. It is shown how this information can be deduced 
fron the case of infinitesimal-amplitude oscillations. The bifurcation theory 
analysis reveals that when the bifurcation parameter, e.g., the angle of 
attack, is increased beyond a critical value at which the aerodynamic damping 
vanishes, new solutions representing finite-amplitude periodic motions 
bifurcate from the previously stable scene;/ motion. The sign of a simple 
criterion, cast in terns of aerodynamic properties, determines whether the 
bifurcating solutions are stable (supercritical) or unstable (subcriticai) . 

For the pitching notion of flat-plate eirfoilr. flying at supersonic/fcyparccnic 
speed and for oscillation of flaps at transonic speed, the bifurcation is 
subscritnal, implying either that exchanges of stability between steady and 
periodic motion are accompanied by hysteresis phenomena, or that poterriallv 
large aperiodic departures from steady motion nay develop. On the other hand, 
for the rolling motion of a slender delta wing m subsonic flight (wing reel.) 
the bifurcation is found to be supercritical. This and the predicted amplitude 
of the bifurcation periodic notion arc m good agreement with experiments. 



£1. UiniODUCTIO.S 




Problem of aerodynamic stability of aircraft flying at small angles 
of attack hove been studied extensively. With increasing angles of attack the 
problem become core complicated and typically involve nonlinear phenomena such 
as coupling between nodes, amplitude and frequency effeers, and hysteres.s. 

The need for investigating stability characteristics at high angles of attack 
was clearly desonstrated by Orlik-Huckenann [1] in his survey paper which 
largely deals with experiments. 

On the theoretical side, the greater part of an extensive body of 
work is based on toe linearized theory, in which the unsteady flow is regarded 
S3 a snail perturbation of sene known steady flow (possibly nonlinear in, 
the angle of attack) that prevails under certain flight conditions. The question 
of the validity and limitations of such a linearized perturbation theory is of 
fundamental importance and yet has been investigated only rarely. One nay ergo* 
that in principle, it should be possible to advance to higher and hir'n'-r 
of attack a by a series of linear perturbations, tince .:he solution ..t each 
step should include a steady-state part which, when added to r.he previo.u 
steadi-state solution, would provide the starting point £< r 
the next perturbation. This may util be true provided that at cj eh step the 
steady motion is stable both statically ana dynamically , and that the actual 
disturbances, e.g., the amplitude of oscillcrirn, remain sccll. However, when 
the angle of attack exceeds a certain critical value at which tr.t ster.Jv 

motion is no longer stable, the linear theory predicts an exponential 6 iovth ot 
the perturbation with tint and, therefore, must itself ce..;e to be v«.Ild .itcor 
finite tine. The notion of the aircraft under these cond'tiono -on only be 
studied using a nonlinear theory. 





jsjoTtjsx: 
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In this paper we investigate the stability charcctorintics of an 
aircraft trimmed to a uaan angle of attack a near a cr ct which the steady 
uotion beccnss unstable. Padfield {2] studied a similar problem, using the 
method of multiple scales, which is valid only for weakly non-linc&r oscilla- 
tions. We shall study the problem by rscanc of bifurcation theory. This will 
allow us to draw on recent mathematical developments (c.g., [3]) that are 
particularly well suited to investigating fundamental questions in linear and 
nonlinear stability theory. A numerical scheme based on bifurcation theory 
was proposed earlier [4] for analyzing aircraft dynamic stability in a rather 
general framework. More recent work by Guichetcnu [5) demonstrates the con- 
siderable potential of bifurcation theory in flight dynamics studies, particu- 
larly toward establishing a uathod for the design of flight control systems to 
ensure protection against loss of control. On the ether hand, while acknowledg- 
ing the importance of the aerodynamic model in determining the aircraft stability 
characteristic!., neither of these works contains an adequate assessment oi Uv. 
model requirements. The treatment or unsteady flow effects, in particular, 
receives no attention. In contrast, we shall focus on just this aspect of the 
problem at the expense of narrowing the scope of the motion analysis. 

We shall restrict ourselves to the simpler case of single-degree-of 
freedom cotions, e.g,, pitching or rolling, of an aircraft about its trimmed 
flight condition. This will enable us to analyze notions for which complete 
aerodynamic information is available, for certain shapes, in the form of exact 
analytical or numerical solutions {6-12] of the unsteady inviscid flow equations, 
or empirically from exrerinents t 13 ] - In this way it will be possible to 
establish a fora revealing a precise analytical relationship between the basic 
aerodynamic coefficients and the characteristics of the motion. 
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Specifically we shall consider the following throe types of oscillatory 

notions: 

A. Pitching supersonic/hypersonic aerofoils in rectilinear flight (Fig. 1); 

B. Flap oscillations in transonic flight (Fig. 2); 

C. Wing rock of slender delta wings in subsonic flight (Fig. 3). 
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$ 2 . MATHEMATICAL FOFKULATION 

Let the aircraft be in level, steady flight up until tine t“0 
when it is perturbed frosa its trim position. During the subsequent notion 
the center of gravity continues to follow a rectilinear path at constant 
velocity . For a single-degree-of-freedon oscillatory notion, the 

equation of notion is 


A 

1*4 - Git; A) 
dt 


( 1 ) 


where the angular displacement 5(c) is the instantaneous pitch angle in 
example A , flap deflection angle in example B and roll angle in example 
C , I is the respective moment of inertia, while G(t) is the corresponding 
instantaneous moment of the aerodynamic forces. 

In Eq. (1), X in general represents a set of parameters defining 
the steady flight at the tram condition, e.g., flight Mach number , ratio 
of specific heats \ , angle of attack a , etc. In this paper we shall 
consider X ro be the angle of attack a in examples A and C , and the 
mean flap deflection angle 6^ in example B . In other words, all the other- 
defining parameters will be held iixad when considering the consequence oi 
varying X on the aircraft motion characteristics. ’Jo assume that tne cement 
required to trim the aircraft at X has been accounted for, so chat G(r;X) 
is a measure of the perturbation moment only. 

• 

It is clear that the instantaneous motion state, C(t) and £(t) , 
and the instantaneous moment G(t) at time t are. a result of the interaction 
of aircraft notion and the unfeady aerodynamic forces from time zero to time t . 





Consequently, the instantaneous moment G(t) depends r.ot only on the 
instantaneous motion state, £(t) and £(t) , but also on the paat notion 
history from time zero to t . This is to say that G(t) is a functional 
of SCtj) , (OSt^St) as shown by Tobak et. al. [ 14 j . Thus 


G(t) - G 5(t 1 ) 


Such a functional, which includes all tine-history effects, is however rarely 
known. 


Now for notions for which £(t^) analytic, the functional is 
equivalent to a function of an infinite set of variables, i.c. 


Gtt) =■ G s(t ]L ) 

V°J 


G(f(t), 5(t), c(t\ 


d n C(t) . 

, - t . . 

uc 


Fot most problems encountered ia the study of the dynamic stability of an 
aircraft, the notion is slow although its amplitude nay be finite or largo. 
Under these conditions S(t) , S(t), ... in Eq. (3) may be neglected and, 
as a first approximation, we get 


G(t> ** G(5(t),£(t)) 



7 . 


Wo further assume that G is an analytic function of 5 and 5 

•2 

Expanding (4) as Taylor series and neglecting terms 0(5 ) and higher for 
slow motions we get, after re-introducing the parameter X , 


5(tH 

G(t;\) - C 0 (Ut);X) + — y- 


(5) 


where !. is a characteristic length and 


C Q (0;X) - 0 


( 6 ) 


as required at the trim condition. We note that the form of the instantaneous 
moment G(t;X) in Eq. (5) is consistent with the exact analytic solution for 
the pitching moment in example A [Ref. 10], the numerical solution for the 
hinge raomenr in example B [Ref. 12] and Che experimentally derived empirical 
formula for the rolling moment in example C [F.ef. 13]. 

Summing up, the singlc-degree-of-f recdom notions considered in tlii., 
paper will be based on the following mathematical problem 


,d Z E 


Ql 


i~ » c 0 (C;X) +— G (5;X) i G(5;5;X) 
dt^ 00 


(7a) 


5(0) 


(710 


5(0) 


(7c) 


The functions G Q (5;X) and G^(5;X) are generally nonlinear in 5 and have to 
be determined from the study of unsteady aerodynamics, either theoretically 
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or expericentally. Evidently Gg is related to the restoring moment and 
G^ to the damping moment. 

In many situations, it is known that when the parameter X reaches 
some critical value X , the aerodynamic damping G^ vanishes and the 
steady flight at the trim condition X loses its stability. The main 
purpose of this paper is to use Hopf bifurcation theory to determine the motion 
characteristics of the aircraft whose trim condition i3 near or beyond X^^ . 
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§3. BIITRCATIOH THEORY 


We introduce the dimensionless time r u V t/Z , where l it a 

co 

characteristic length, equal to the chord length o£ the aerofoil in example 
A, the chord length of the ilap in example B and the chord length, of the wing 
in example C. Hereafter we use (*) to denote d/dr . We furtner let 


F(C ;£;>.) - G(C,?;X)/(I vjv* 2 ) 


F 0 (5;*> + £F^(C;X) 


then Eq. (7a) may be written 


li - I 

dt ^ 


- F (?,e;A) 


An expansion of F(5,5;X) in a Taylor series in f end C and 
change of notation ^ « £, u 2 “ ? yield for Eq. (9) 


h " A..(X)u. + E ijk (X)u jUk + C.. kt (X)u jVjl 


+ 0( | u| "*) , (i - 1,2) 


where 


0 1 


-S(X) -D(X) 


S(\) " -F 0 '(0;X) , B(X) « -F^OjX) 


3 ljk “ ° ’ B 2jk " 2! 3 u.3u"L 

2 u“0 
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(Although Eq. (9) has been derived on the assumption of slow oscillations, 
our subsequent bifurcation analysis of (9) will hold for general f(€,(;X) , 
i.e., as if no restriction had been placed on the magnitude of l .) 

In Eq. (11) the tensors 3 and C represent the effects of finite 
amplitude to the second and third order. We note that the following symmetry 
properties hold: 


D 2jk “ B 2kj 


C m C ■ C . m c o c 
2jk l 2 j 5.1; u 2kxj 2£jfc 2lkj 


(12a) 

(12h) 


On the basis of Eq. (10), ve shall 3tudy now Ihe linear and nonlinear 
stability of the motion. 

3.1 LINEAR STABILITY THEORY. 

The stability of the steady motion at the trim condition X to 
infinitesimal disturbances is determined by the nature of the eigenvalues of 
the matrix A. They are 


V X) " I t~ D(X > 1 *4> 2 (X) - 4SU) 1 (13) 

Case I : S(X) < 0 . In this case n^ > 0, n ? < 0 . The steady 
motion at this condition X is '-.lways unstable. 


• -T-- 
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Case II ; S(X) > 0 . 

Ila : D(X) < 0 . In this case Re(n^) > 0 and the steady 

motion at X is unstable. 

lib: D(X) > 0 . In this case Re(n ) < 0 and the steady 

1 

2 

motion at X is stable. 

Thus, only in Case lib, when both stiffness derivative S(X) and 
damping derivative D(X) arc positive, is the steady motion at the trim 
coididon X stable to infinitesimal disturbances. In fact, stability theory 
13] can be used to show that stability of the steady motion in this case is 
assured if and only if the disturbance is sufficiently small. 

In all cases in which tha linear theory predicts growth of the 

disturbance amplitude, the growth predicted is of exponential form and 

hence the linear theory au'-.t cease to be valid after come finite tine .dun 

the amplitude is no longer small. Thus, what eventually happens to a not r or. 

for which linear stability theory predicts a growth of disturbances cannot 

be determined from the linear theory itself. Instead, the full nonlinear 

inertial equations of morion, or a suitable approximation of them, such as 

(10), must be adopted to determine tne ultimate state of the motion. Of 

particular interest is the dynamic stability boundary X « X , whert 

S(\ X > 0 and, D(X ) = 0 . The stability characteristics near this 

cr cr 

boundary will be studied presently. 

3.2 HOFF BIFURCATION THEORY 

At the dynamic stability boundary X » X^ , u. have S(X^ r ) > 0 



( 1 / 


n (x 

i 


cr 


asslx 7 = ±io 

cr 


0 


The existence of purely imaginary eigenvalues of the m«trix A at X « X cr 

is the characteristic sign of a Kopf bifurcation [3,15], signaling a changeover 

from stable steady motion to periodic motion. On crossing X •» 1 cr , the steady 

motion that had been stable for X < X will become unstable to disturbances, 

cr 

resulting (after a transient notion has died away) in the existence of a new 
motion, vhicn (if it is stable) will be periodic. In the vicinitv of X“X , 
the circular frequency of the periodic motion will be nearly equal to . 

We coll the new solution of the equations of motion a bifurcation solution. 

In this section we shall determine its character and a criterion for its 
stability. 


For X slightly larger tnan X , the eigenvalues of the rntrix A 

are 

n - - i D(x) t in(x) (15) 

l 2 

n 

where 

Q(X) - /s<X) - D 2 (X)/4 (16) 

We shall assume that 

D'(X } < 0 (17) 

cr 

which is the usual case in applications. (The case D*(X ) " 0 can be 
created in exactly the uame way.) The normalized eigenvector 5 (X) 



I 



t 1 * ' 


13 . 


associated with the eigenvalue n(X) is 


C(X) 


C j(X) 

1 - i 


C 2 ( x > 

I < 


n(X) 


(18) 


whereas the adjoint eigenvector Q"(\) with eigenvalue n(X) , which is 
the cosplcx conjugate of n(X) , is 


C*(X) 


* ’ 
;?(x) 

1 +i 

’ » 
n(X) + D(X) 

C$(X) 

2nKx7 

1 


(19) 


A. Hopf Bifurcation. The bifurcation solution u(t t X) aay bo 


written as 


u « a(r); + a(r)c 

Following Iooss and Joseph ([3J, p. 125), we get 


CO) 


a •* eb^s) + e^'o 2 (s) + c^b^is) t 0(c ’) 
2 ^ 

s ■ [u Q + + 0(c ) jt 

X « X + c^X. + O(c^) 
cr 2 . 


( 21 ) 


where, for brevity, we cait the lengthy solution foras .'or r, b , u^, , 

and X 2 (cf. (HD* The solution is periodic in t vi h circular frequency 
2 4, 

equal to Uq t t o, + 0(c ). 
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B. Stability o£ the Bifurcation Periodic Solution. According to 
Fioquet theory [3], the stability of the bifurcation periodic solution (21) 
is determined by the sign of an index u . To 0(c ), u has the fora 

u - D'(X cr n 2 c 2 4- OCrS (22) 


and the bifurcation periodic solution is stable if u < 0 , unstable if 
n > 0 . Since we have assumed D’(\ cr ) < 0 » stability thus depends on the 
sign of , with X^ ■» 0 denoting stability and X^ < 0 instability. It 
remains to cast X^ in raore recognizable terms. After considerable 
manipulation, we get 


•> 

c“ 


f 

2 2 
3 r ^ 2 3 F 

3 2 F 2 

3 J F 2 S't! 


■ , 3 
4w 0 


- 2 + “o . 2 

3u^ 3u, 

3u^3u£ U 0 

, 2- ^ u 0 . 3 

[ 3u 1 3u 2 3u 2j 

l 


(23) 


in terns of F(u^,u,;X) . Fron (8) we sec that the function ? is directly 
related to the moment C'5,F,;X) acting on the aircraft which is performing 
a f initc-acplitude oscillation C around the trim condition X . Equation 
(23) deconstratos that the stability of the periodic motion near the dynamic 
stability boundary X is determined by the behavior c£ the aerodynamic 
response G(^,C;X) in that vicinity. 


With the assumption of slow oscillations under which the fora of 
(5) was derived (terns of 0(5 2 , O neglected), we nay substitute (5) with 
(8) into (23) to got 


u 



(T "r * 
U 0 * 1 


2„ , 

u 0 ‘l 


F,“0 


cr 


(24) 
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Since froa (14) and (lib) 


“0 " ’ • C *7 V! '\r > 


Eq. (24) becomes, after using (8) 


or, 


o 

3 

N 

U 

A . . 


[1 

4 

dS 

. F o' K ' X =r>. 

1 


e-o 


u " - 


2 

c u o 

d 

' V 

<^ ) cr )1 

“ 

4 

dC ' 

.V 

« = X cr>. 

_1 


Using (10) we get 


C 2112 ( V,"' 

2 s ‘ x ,r' 


We have thi s established the following criterion: the bifurcation perioi 

notion is stable or unstable according to 

y < 0 or u > 0 


or alternatively, since 5(1^) > 0 , according to 


C 2112 <X cr ) < ° or ^ll^cr 5 > 0 


(25) 


(26 a) 


(26b) 


(26-) 


(27a) 


(27b) 


The two possibilities are well illustrated in the fora of bifuicaticn 
diagrams as 3'novn in Figs. 4a and 4t>. 
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In a bifurcation diagram, the abscissa is the bifurcation parameter 

X , while the ordinate is a parameter characteristic of the bifurcation solution 

alone. In our casa it is t , a noasure of the amplitude of the bifurcation 

periodic solution. Stable solutions are indicated by solid lines, unstable 

solutions by dashed lines. Thus over the range of the bifurcation parameter 

X < X where the steady-state notion is stable, e is aero, and the stable 
cr 

3tcadv notion it represented along che abscissa by a solid line. The steady 

notion becomes unstable for all values of X > X as the dashed I’ne along the 

abscissa indicates. Periodic solutions bifurcate froa > n X either super-* 

cr 1 

critically or subcricically. 

When C_, (X ) < 0 hence u < 0 (implying \ n > 0) , the 
bifurcation is called supercritical and its characteristic fora is shown in 
rif. 4a. In this case, stable periodic solution-, (solid curves in Fig. 4o) 
exist for values of X > X . The amplitude of the periodic solution at c 

given value of X - X is proportional to c , hence is vanishingly small unei 

1 

X ~ X is small, varying essentially as (X - X c _) . 

When > 0 hence t > 0 (implying *■ O'* , the 

bifurcation is called subcritical and its characteristic form is shewn in 

Fig. 4o. In this case, periodic solutions exist for values of X < X , 

but they are unscaole (dashed curve in Fig. ^b). Whether 3toble periodic 

solutions do or do not exist for X > X depends predomir antly on the 

behavior of the danpirg G^(C;X) for X > X . If r.o sura stable periodic 

solutions exist for X > X , then when \ is increased bevond X 

cr - cr 

aircraft may undergo an aperiodic motion whose departur? fisc the steady 




the 
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notion nt 1 * X is potentially large. 

In the core likely event that stable periodic solutions do exist for 

X > X , their amplitudes must be finite, and net infinitesimally small, even 

for small positive values of X - . It is likely that this branch of stable 

periodic solutions will join that of the unstable branch in the way illustrated 

in Fig. 4b. In this event, the fora of the bifurcation curve for values of 

X < \ cr helps explain the situation where the steady-state motion could bp stable 

to sufficiently small disturbances but become unstable to larger disturbances. 

Thus, Fig. 4b suggests that for X < X , so long as disturbances arc of small 

cr 

enough amplitude to lie veil below those of the unstable branch of periodic solutions 

(curve OB in Fig. 4b), they will die out with time and the steady motion vs! ) 

remain stable. However, disturbances with amplitudes, sufficiently larger t»u;r> 

those of the unstable branch nay actually grow up to the ultimate notion .,0 

t -* « , which will be that of the stable blanch of periodic solutions (cu.vt 

BA in Fig. 4b). Fins’ ly, we note that if the motion does attain the 

branch of periodic solutions (say, for X < X^ ) then hysteresis effects will 

manifest themselves with further changes in X . When X is increased 

beyond X , the motion will continue to be periodic with finite amplitude 

(point A in Fig. 4b). If X is now decreased below X , the periodic 

motion will persist, even at values of X where previously there had been 

steady motion when a was being increased. Hot until X is diminished 

beyond a ceitriin point (point B in Fig. lb) will the notion return to the 

steady- state condition (point C in Fig. lb) that had beer, experienced when 

\ wa3 increasing. 
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54. APPLICATIONS TO AIRCRAFT NOTIONS 


In this section we shall apply the thtoiy developed in 55 2 and 3 
to three different sinele-dence-of-freedoa motions of the aircraft. 

4,1. PITCHING MOTION OF SUPERSONIC/HYPERSONIC AIRFOILS 

In the case of pitching oscillation (Fig. 1) of a supersonic/ 
hypersonic aerofoil in rectilinear flight, £(t) in Eq. (1) is the instantaneous 
pitch angle G ( t) measured from its trim condition, \ is the mean angle of 
attack a , and G is the pitching moment about the center of gravity. Thus 


G(t;X) •* qSt C (t; a, M^, Yt h) 


( 18 ) 


where q iG the dynamic pressure, and S and Z are the reference arm 
and length. The instantaneous pitching moment coefficient C^ is in general 
also a function of the angle of attack, the flight Kr.ch number M , the 
ratio of specific heats y , and the pivot axis position h . 

For large amplitude slow oscillations exact analytic solutions of 
C^ exist for simple shapes [10,11] and they arc in the fom consistent with 
Eq. (i), i.e. 


C (t;«) - C (e(-0;u) t 6 C (0(r); a) - C (0;a) 


m 


(29) 




Moreover, it is shown that [10,11] 


C (G;a) - C (0-ra) , (i«-0,l) 

m. n. 


(30) 


••Cm. 

i 


Accordingly 


30 


3Cu. 

i 

3a~ 


(i“0,l) 


(31) 



and Eq. (26b) reduces to 


y 




(32) 


where the stiffners derivative S(a) and the dauping-in-pitch derivative 
D(a) ct angle of attack a are related to the pitching moment coefficient 

c C! b y 


and 


S(cs) ■ - C' (a) , D(o) » - C (a) 

D>- u. 


D(a ) *» 0 
cr 


(33) 


(34) 


A typical example of the damping derivative D(a) versus angle of attack a 
is shown in Fig. 5 [16]. 


Wo therefore conclude that when the angle of attack u is increased 
beyond a critical v.ilte a ^ at which the aerodynamic damping vanishes, i„c., 
D(a cr ) ■ 0 , Che steady flight at a cj _ loses its stability and, after a 
transition, results in a iinite amplitude periodic notion. Furthermore, the 
stability of this bifurcation periodic motion depends on whether D'(u)/S(ct) 
increases or decreases on crossing the stability boundary . 

By utilising an approximate relation [17] 


D(a) « b[S(a cj .) - S (a) ] , b > C 


in (32), we get 


2 , 
e u Q o 


— ^ Ski S(o) j 
da | a «■ a_ 


cr 



For the cast cf a flat plate aerofoil in supcrscnic/hypersonic flow 
the stiffness and damping-in-pilch derivative, S(e) and D(a) , arc known 
exactly in analytical form [7] for the angle of attack up to the shock detachment 
angle. It is shown in Table 1 that u > 0 for all combinations, of M w and h , 
so the bifurcation is subcntical and the bifurcation periodic notion ic unstable, 
implying that exchange of stability between steady and periodic motion are 
accompanied by hysteresis phenomena, or that potentially large aperiodic 
departures from steady flight may develop. 
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4.2 FLAP OSCILLATION IN TRANSONIC FLOW 

The case of a flap oscillating about the hirge is conceptually similar 
to that of an aeroloil pitching about a pivot axis. Thus, the oscillating flap 
may be regarded as an aerofoil pitching in a non-uniform incoming stream that 
results from a uniform free stream passing the fore-body fixed aerofoil. (Fig. 2) 

In the present case £(t) in Eq. (1) is the instantaneous flap 

deflection angle 6(t) , X is the mean deflection angle 6 , and G is the 

m 

hinge moment. Thus 


G(t;X) - q Si (ytjd^.M^.Y) (3S) 

where is the instantaneous hinge moment coefficient. In the uuLhet atical 
modeling of [12], the Cj used is of the form 


C. (i;6 J » C. (6(t); 5J + 6 C, <6( r);6 ) - C (0;6 ) (36) 

h m 11 q m h^ tn m 


Eq. (36) is consistent with (5) and was validated in [12] by comparisons with 
results of large sea’c numerical integration of the coupled inertial/f icwfield 
equations. 


It should be pointed out that the form (5) or (36) for the moment 

coefficient enable, one to calculate C, (6 (t); 6 ) by linearising, for small 

amplitude ar.d frequency, the flowficld equations about the steady flow 

corresponding to 6 and 6(t) . Such calculations of C, (6(t);6 ) from 
m m 

the linearised equations require negligible amount of computing time compared 
with the nonlrncar method used in [12j, yet, they are consistent with the level 
of approximation leading to (36) or (5). With this improvement on the method 
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cE calculating C. (5(t);6 ) , the main conclusion of f 12 ] is further 

ns 

strengthened in that it costs very little, in all cases, to calculate the notion 
of the flap using the mathematical modelling approach of Tobak at. a. [ 14 ] than 
to solve the coupled iner tial/f lowf ield equations. 

It was further established [11] that 


C, (6 j 6 ) « C. (6 + 5 ) (i*=0,l) 

n. m n. n 

i l 


(37) 


which in turn reduces Eq. (26b) to 

2 r 


e u„ 


, D*(« ) 

d u 


d6 SE6 ) 
m n 


6 “6 


- 1 m cr 


where 


(38) 


S(c ) *= -C, ' (5 ) , D(6 ) - -C. (o ) , D(6 > - 0 

m hp a n n cr 


(i9> 


We therefore reach a similar conclusion regarding the notion characteristics 
of the flap oscillation in the neighborhood of 6 35 t’ at in i4.1 regarai \g 

the notion characteristics of the aerofoil pitching motion in the neighborhood 
of a 

cr 


Numerical results of Ref. 12 for S{6 ) end D(6 ) are reproduced, 

m m 

using spline fitting, in Fig. 6 tor N/.CA 64A010 aerofoil at aero mean angle of 
attack in a transonic stream with M *■* 0.8 and y » l.a, It is found tlat 


dfi 3 (.6 ) 
in G 


= -59.20 end hence the bifurcation is subcritical and the 


6 *=6 


- 1 m cr 


bifurcation periodic solution is therefore unstable. Th-'s implies that 


similar 
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to the flat plate at superporic/hypersonic flow, exhange of stability betvreen 
steady and periodic notion o£ the flap of the transonic aerofoil are accompanied 
by hysteresis phenomena, or that potentially large aperiodic departures from the 
strady potion may develop. 

4.2. WING ROCK OF SLENDER DELTA WINGS II! SUBSOHIC FLIGHT 

The phenomenon known as "wing rock" of a slender delta wing in 
subsonic flight has been a subject for intensive investigations by many researchers 
[13, 18-21) and is now well documented. It is known that the steady flow past 
a slender delta wing at small enough angle of attack a is syssaetric and 
remains steady. However, as the angle of attack a is increased past a critical 
value the nyrmetric configuration of the lending edge vortices becomes 

asymsatric, causing a loss of roll damping at small angles of roll [21]. 
Consequently, at a S a ^ small disturbances introduced into the flowfield cause 
the rolling motion of the wing to develop, resulting ir. wing rock. 

For rolling motion of the wing, £(t) in Eq. (1) is the instantaneous 
roll angle $(t) , 1 the steady angle of attack a , and G the rolling moment. 
Thus 

G(t;\) “ qS£, C ? (t:a,A) (40) 

where, for incompressible flow, the instantaneous rolling resent coefficient 
depends on the angle of attack a and the sweep back angle A of the 
delta wing. It is .ound from experiments [13,18] and numerical calculations 
[19] that for a given A a good approximation to the instantaneous rolling 
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moment C f is of the form 


C. » C. («(t);o) h $ C U(T);a) 
1 *0 1 


( 41 ) 


which is consistent with Eq. (5). 


An an example, for the case of an 80” sweep-back flat delta wing 
ef. 18, is well approximated by a 

to the following equation of the rolling motion. 


used in Ref. 18, C is well approximated by a power series [19] which leads 


<? “ F(4>,<p;a) 

- [b 1 Cct)<J> + b 3 (e)* 3 ] + «[b 0 + b 2 (a) b 4 (a)4> 2 ] 

S F Q (();a) + Fj^jcO 

where, with different scalings in [19] accounted for. 


(42) 


b. (a) ■ < C. a. (a) 

i 1 J- 


bj (a) « x C x a^-(a) 


(i-1,3) 

(j-2,4) 


(43a) 


(** 


43b) 


In (43), x is a factor arising from the different scalings used for the time 

variable in Ref. 19 and in the present paper. Thus, from [19] 
x - 2c /L - 2 x 0.429/0.107 4 8, ^ - 0.CJ8 and a. (a) are tabulated in 

Table 5 of Ref. 19 which yield the following table for b.(ct) 
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a 

10° 

15° 

20° 

25° 

b i 

-0.0265 

-0.0721 

-0.1S77 

-0.3320 

b 2 

-0.0101 

0.0090 

0.0596 

0.0959 

b 3 

-0.1222 

-0.2714 

-0.0501 

0.2894 

b 4 

0.1491 

0.1159 

-0.1799 

-0.9977 


From Eq. (lib), the linear damping in roll is 


D(a) " -F^^a) “ ~b Q - l" 2 (a) 


(44) 



i 


t 


& 

t 

P 

ii 



j 

n 

j 

x 


} 


where -bg >0 is proportional to Lite damping coefficient for the bearings 
in sting :n the experiments of [18]. In order to compare our theoretical 
prediction with the experiments o £ [18] we choose the value of -bg ruch 
that Eq. (44) yields zero damping in roll at the same angle u 4 18.6° 
as that in [18]. The function D(a) is then plotted in figure 7 using splint- 
fitting. 


From Eq. (26c) the index for stability of the bifurcation periodic 
rolling notion is 


c b 4 <c tr ) 

2 bjTTT 


(45) 


At a “ 18.6° , we know from the above table that 
cr 

also find from the spline-fit curve for b^(a) (Fig. 


b, (a ) « ~ 0.1591 
l v r 

8) that b. (a ) “ 
4 cr 


0 , we 
0.05473 < 0 


hence p < 0 , and the bifurcation is supercritical, implying that the 




bifurcation periodic rolling notion i3 stable, in agreement with [18]. 


We now further compare the supercritical bifurcation diagram 
predicted by the present theory with experimental results of Ret. 18. Combining 
(21c) with (22) to eliminate X, and then use Eq. (45) tc get 


a - a 


cr 


- c 


h. 

0 

2 D'b, 


cr 



(46) 


The last equation was obtained by using Eq. (25) and (42) . From Fig. 7 i-ttl 
noting that 1° ™ u/180 , we get D' ( Q cr ) = -0.6131. Hence Eq. (45) re*' -c< ; 
to 

a -a =■ C. 1120 e 2 (47) 
cr 


where c is the amplitude of the bifurcation periodic solution [11], 

Eq. (47) is plotted in Fig. 9 for the 80° sweep back delta wing and compared 
with experimental results of Levin and Katz [18] . The agreement is seen to 
be excellent in the neighborhood of where the bifurcation theory applies. 

It should be noted that although u < 0 and hence the bifurcation 
is supercritical, the /alue of |vj is snail in the example. This incites 
that the bifurcation is close to the boundary between subciitical and super- 
critical. Consequently, on .ncreasing the angle of attack a by a small 
amount the amplitude o'. tne resulting periodic notion will be quite large. 
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This theoretical prediction is confirmed by experiments as seen from Frg. 9; 
for instance, increasing the angle of attack by 1° past « cr (-18.6°) results 
in a * table periodic motion whose amplitude is 22.6°. 





55. CONCLUSIONS AND DISCUSSIONS 


We have shown how bifurcation theory can be used to study the 
nonlinear dyrsmic stability characteristics of an aircraft subject to single- 
degree-of freedom motion about its trimed steady flight near the stability 
boundary. In most cases of slow notion, the required moment of the aerodynamic 
forcer, in in the fora of Eq. (5). The theory shows that when the bifurcation 
parameter A , c.g«, the angle of attack, is increased past the stability 
boundary X where the aerodynamic damping vanishes, the steady flight loses 
its stability, rcsulting(after the transient motion has died out) in a finite- 
amplitude periodic motion. We have also established a simple criterion (Eq. (27)), 
for the stability of the bifurcation periodic motion in terms of the aerodynamic, 
coefficients. The theory predicts that the bifurcation solutions are unstable 
(subcritical) in examples A and 3 and stable (supercritical) in example C, tho 
latter prediction being also in good agreement with available experimental data. 

In the case the theory predicts subcritical bifurcation, the abrupt 

change resulting from increasing X past X could cause an abrupt structural 

cr 

change of the flow fieM which may in turn r.- ider invalid the font of the 
moment of the aerodynamic forces, Eq. (5). Under tense conditions cercdyriamic 
information in a different form cay be needed. However, the theory as developed 
ir. this paper is valid up to A and can be used to predict the onset of 
subcritical bifurcation. 

In the case the theory predicts superc -itical bifurcation, the 
bifurcation periodic solution is stable to small enough disturbances for \ 
in the neighborhood ol A . However, with further increase in X , the 
periodic nocion might lose its stability, causing mother bifurcation at 
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X * > X^ r , oay, which can be either subcritical or supercritical. The 

resulting bifurcation solution nay also be almost periodic and this sequence 
of bifurcation nay continue. 

Tinnlly, the theory in principle nay bo generalized to the notions 
of the aircraft involving more than one-dcgrec-of-frcedo-n by the method of 
projection (3). In that case one should be aware of the possibility of chaotic 
motion (strange attractor) occuring after a finite number of successive 
bifurcations. 
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Table 1 Values nf Gtability criterion V»(H OT , h) for flat-plate airfoil: 
c ™ 1, Y » 1.4; w > 0 oubcritical bifurcation, p < 0 auporcritical 
bifurcation 


0.1 0.2 0.3 0.4 



0.1 0 .? 0.3 0.4 



* 1 
w 

? 

1.5 

36.5 

19.0 

-- 


— 

2 

39.5 

25.8 

17.2 

13.0 

14 . 7 



1.6 

39.6 

23.4 

13.6 

8.8 

10.7 

3 

57.2 

38.6 

26.8 

20.9 

23.3 

./ 

'1 

» * 
.< * 

1.7 

40.0 

24.8 

15.5 

11.0 

12.8 

4 

84.7 

58.4 

41.4 

32.8 

36.4 

' 


1.8 

39.7 

25.3 

16.3 

12.0 

13.7 

5 

107.3 

74.8 

53.7 

42.9 

47.5 


r , i t 

S ,' 

1.9 

39.4 

25.5 

16.8 

12.6 

14.3 

6 

123.6 

86.7 

62.6 

50.3 

55.6 


2.0 

39.5 

25.8 

17.2 

13.0 

14.7 








-* 

\ h 













") 


0.25 

0.26 

0.27 

0.28 

0.29 

0.30 

0.31 

0.32 

0.33 

0.34 

0.35 


1 














k ” 

l ; 

l 

1.6 

] 0.6 

10.1 

9.7 

9.3 

9.0 

8.8 

8.6 

8.5 

8.4 

8.4 

6.5 

- 


1.7 

12.6 

12.2 

11.8 

11.5 

11.2 

11.0 

10.8 

10.7 

10.6 

10.6 

10.7 


/ 

/ 

1.3 

13.6 

13.2 

12.8 

12.5 

12.2 

12.0 

11.6 

11.7 

11.6 

11.7 

11.7 



1.9 

14.1 

. 13.7 

13.4 

13.0 

12.8 

12.5 

12.4 

12.3 

12.2 

12.2 

12.3 

! 

I 

r 

l 

?.o 

14.5 

14.1 

13.8 

13.5 

13.2 

13.0 

12.3 

12.7 

12.6 

12.7 

12.7 

x~vr v*.<a WM- + 



31 


ACKNOWLEDGMENTS 

This research was financially supported by a NASA Grant NAGW-575 
for which W.H.H. is grateful. 


REFERENCES 

1. K.J. Crlik-Ruckcmann, Dynamic Stability Testing of Aircraft - Needs Versus 
Capabilities , Progress in Aerospace Sciences, Vol. 16, Ho. A, Pergamcn 
Press, New York, pp, 431-447 (1975). 

2. G.D. Padfield, Nonlinear Oscillations at High Incidence , AGARD CP-235, 
Dynamic Stability Parameters, Paper No. 31, Hay 1978. 

3. G. looss and D.D, Joseph, Elementary Stability and Difurcation Theory , 
Springcr-Varlag, Hew York, 1980. 

4. 3.K. Kehra and J.V. Carroll, Bifurcation An al ysis of Aircraft lliph An gU - 
o C-Atsack Pligh t D ynaai cs, Nov Approaches to N onlinear Problems in 
Dynamics , Ed. P.J. Holmes, SIAM, Priadelphia, pp. 127-146 (1980). 

5. P. Guicheteau, Bifurcation Theory Applied to the Study of Control Loss es 
cn Conaat Aircraft' , La Recherche Aerospatiale, pp, 1-14 (1982-3). 

6. H.H. Hui, Stability of Oscillating Hedges and Caret Wings ir. Hyp ersonic 
and Su p et sonic Flow s, AIAA J., Vol. 7, No. 8, pp. 1524-1530 (Aug. 1969). 

7. W.H. Hui, Supersonic/Hypersonic Flow Past an Oscillating Fl at-Plate at 
High Angles of Attack , ZAMP, Vol. 29, Tasc. 3, pp. 416-427 (1978). 

8. H.H. Hui, An Analytic Theory of Supersonic/Hypersonic Stability at High 
Angles of Attack , AGARD CP-235, D ynamic Stability Parameters , Paper No. 22, 
Kay 1978. 

9. W.H. Hui and M. Tobak, Unsteady Newton-Eusenann Flow Theory . Part I: 
Airfoils , AIAA J., Vol. 19, No. 3, pp. 311-318 (Kar. 1931). 



32 



(}? 



:<l * 


i 


l ' 1 





;* 




10. W.H. Hui, Large-Amplitude Slow Oscillation of Wedges in Inviscid Hyp e rsonic 
and Supersonic Flows , AIAA J., Vol, 8, No. 8, pp. 1530-1532 (Aug. 1970). 

11. W.H. Hui and M. Tobak, Bifurcation Analysis of Aircraft Pitching Mot-ions 
About large Mean Angles of Attack . J. Guidance, Control, and Dynamics, 

Vol. 7, pp. 113-122, 1984. 

12. W.J. Chyu and L.B. Schiff, Nonlinear Aerodynamic Modeling of flap 

Oscillations in Transonic Flow: A Numerical Validation, AIAA J., Vol. 21, 

No. 1, pp. 106-113 (Jan. 1983). 

13. L.T. Nguyen, L. Yip, and J.R. Chambers, "Self-Induced Wing Rock of Slender 
Delta Uing3," AIAA Paper 81-1883, Atmospheric Flight Mechanics Conference, 
Albuquerque, KM, Aug. 1981, 

14. K. Tob&k and L.B. Schiff, The Role of Tiee-Hi story Effects i n the 
Formulation cf the Aerodynamics of Aircraft Dynamics , AGAFD CF-235, 

Dynamic Stability Parameter s, Paper No. 26, May 1978. 

15. E. Hopf, Abauaigung einer periodischen Losung ven einer Statiorf i ren lie- eng 
einea Differentialsystena , Berichten der Mathematisch-Physischen Klaase 
der ~achsischen AkaJemie der Wissenschaften su Leipzig, Vol. XCIV, pp. l-?2 
(1942). 

16. W.H. Hui, "Unified Unsteady Supersonic/Hypersonic Theory of Flow pa a t 
Double Wedge Aerofoils", Journal of Applied Mathematics and Physics 
(ZAMP) vol. 34, pp. 458-488, (1983). 

17. M. Tobak and L.B. Schiff, On the Formulation of the Aerodynamic 
Characteristics in Aircraft Dynamics , NASA TR R-456, Jan. 1976. 

18. D. Levin and J. Kata, "Dynamic Load Measurements uitn Delta Wings 
Undergoing Self-Induced Roll-Oscillations," Journal of Aircraft, Vo 1 . 71. 
pp. 30-36. (1984). 


/ 



I 




1 

K* 


19. K. Constsdinopoulos, D.T. Kook, end A.H. KoyCol., "Sobsonio Vint Kock ° f 

Slender MU «n 8 \ A1AA 85-0198, •)»»»«? 1985. 

20. C:n. »«n end C.E. Un, “theory o£ Kins Reek". AIAA paper 85-0199, 

January 1985. 

21. L.E. Ericsson. "18= Fluid Mechanics of Slender «in £ Rock", Journ»!_£f 
Aircraft, Vcl. 21, up. 322-328, (1984). 


■ n 




\ 



I 


/ 


FIGURE CAPTIONS 
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Fie- 1. Pitching aerofoil in supersonic/hypersonic flight 

Fig. 2. Flap oscillation in transonic flight 

Fig. 3, Wing rock of slender delta wing in subsonic flight. 

Fig. 4. Typical forms of bifurcation diagrams near the dynamic stability 

boundary X where D(X ) =» 0 
J cr cr 

(a) supercritical, u < 0 . 

(b) subcritical, y > 0 . 

Fig. 5. Damping-in-pitch derivative D versus angle of attack a for a 
bi-convex circular arc aerofoil. (Ref. 16) h •* 0.5, y n 1.4 , 
x - 0.075 . 

Fig. 6a. Stiffness derivative S versus mean deflection angle 6 of a 

m 

transonic flap on NASA 64A010 aerofoil (Ref. 12). M^O.S „ y n l.4 , 
a K 0 . 

Fig. 6b. Damping derivative D versus moan deflection angle 5^ of a 

transonic flap on NASA 64A010 a< rofoil (Ref. 12). M n “0.8 , y» 1.4 , 
cs“0 ,6 “2 7 ° . 

Fig. 7. Damping in roll D of a 00“ sweep back flat delta wing versus 

angle of attack a in incompressible flow (Refs. 10,19), Eq. (44). 

a ■■18.6' 1 . 
cr 

Fig. 8. Rolling aerodynamic coefficient b^(ct) of a 80® sweep back flat 
delta wing versus angle of attack a in inccmprcssiole flow 
(Refs. 18,19). Observe b^(o ) ** -0.054/3 «• 0 . 

Fig. 9. Amplitude - of bifurcation periodic motion of a 80" sweep back 
delta ving versus ingle of attack a . 
experiments (Ref. 18). o.^_-ld.6 6 . 


Eq. (47); O 
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